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1 Introduction 

The prediction of risk for ventricular arrhythmias 
after myocardial infarction (MI) remains a clinical 
challenge. Despite various approaches on the basis 
of different invasive and non invasive techniques the 
predictive value of each of these methods is still 
limited. The pathophysiological background reflects 
the combination of an arrthythmological substrate, 
the myocardial scar, and various trigger mechanisms 
like transient ischemia and changes of the activity of 
the autonomic nervous system. The complexity of 
these underlying mechanisms explains the diffi¬ 
culties in the prediction of arrhythmogenic events. 
The quantity of standard procedures such as ECG 
and Holter ECG, late potentials, heart rate 
variability, baroreflex sensitivity and electrophysio- 
logical studies elucidates the confusion and lack of 
sufficient methods in clinical practice [1], Both 
myocardial scar and trigger mechanisms lead to 
inhomogeneity of ventricular de- and repolarization. 
Therefore there is a rationale for exploring new 
avenues such as MCG in the detection of arrhythmo¬ 
genic properties. Various studies have demonstrated 
potential benefit of MCG compared to ECG in 
analyzing cardiac activity, especially when multi¬ 
channel systems with over 60 cha nn els are used. In 
the field of risk analysis various methods of regis¬ 
tration and data analysis seem to be promising. 

In the following, the results of several studies deal¬ 
ing with late ventricular activity, QRS-fragmenta- 
tion, mapping of QRS integrals, relative smoothness 
analysis and QT dispersion will be summarized. 

2 Methods and results 

2.1 Late ventricular magnetic activity 

In MI patients ventricular tachycardias (VT) are 
based on re-entry mechanism s which arise from 
delayed and fragmented depolarization in the 
damaged myocardium. These microsignals, called 
late potentials, can be registered on the basis of 
signal-averaged ECG [2], One of the most important 
contributions leading to widespread acceptance and 
use of this method came from Michael Simson, who 


developed standardized criteria [3], In analogy to 
ECG, Erne et al. first reported on fragmented low 
amplitude magnetic signals at the end of the QRS 
complex based on high resolution MCG recordings 
[4], Two approaches for further analysis were used: 
time domain (TDA) and spectrotemporal analysis 
(STA). 

For TDA a magnetic power spectrum in each 
measurement location was formed from one filtered 
component of the magnetic field using the Hilbert 
transform. Analogous to ECG, the duration of 
filtered QRS complex and of low amplitude signal 
< 1 pT as well as root mean square amplitude 
(RMS) of the last 50 ms of the QRS complex are 
used in order to define the presence of magnetic late 
fields. There is no general conformity concerning 
the cut off points for each parameter. An overview 
of some of the studies dealing with time domain 
analysis is given in Table 1. 


Table 1: Overview of several studies dealing with 
TDA. Sens.= sensitivity: Spec.=specificity; other 
abbreviations, see text. 



Parameter 

Sens 

% 

Spec 

% 

Stroink et al. [5], 
MI,n=12; VT,n=15 

QRSRMS-40<300fT 

67 

67 

Makijarvi et al. [6] 
MI,n=10; VT,n=10 

QRSdur> 115ms 

80 

80 

QRSrmS- 60<1 .200fT 

60 

80 

Korhonen et al. [7] 
MI,n=10; VT,n=10 

QRSdur> 110ms an d 
QRSRMS-40<750fT or 
LAS (i pT)>40ms 

90 

100 

Makijarvi et al. [8] 
MI,n=10; VT,n=10 

QRSdur> 120ms 

70 

80 

QRSrmS-50<7001T 

60 

80 

Korhonen et al. [9] 
MI,n=61; VT,n=42 

QRSdur> 120ms 

55 

90 

QRS d ur> 120ms 
±QRSRMS-40<600fT 

67 

74 




The results for MCG are comparable to signal 
averaged ECG. In the latter, sensitivity varies 
between 63 and 93 % and specificity between 55 
and 76 % [10]. 

Following the idea that spatial aspects of the 
distribution of late magnetic activity might influence 
the results, Montonen et al. compared MCG 
measured from the frontal plane with a measurement 
from the left side of a supine patient and could show 
that the latter improves the identification of patients 
with VT especially in terms of QRS duration [11]. 
The analysis of spectral variability during the QRS 
complex corresponds to the detailed description 
given by Kelen [12] for HR-ECG and can be 
described briefly as the following: First the signal is 
differentiated, comparable to high-pass filtering. 
Second, a sliding-window fast Fourier transform is 
performed over the differentiated signal resulting in 
a spectrocardiogram. Third, the onset and offset of 
the QRS are determined as time instants where the 
total power of the spectrum exceeds a threshold 
value obtained from a noise interval. Finally, the 
variation of the spectra as function of time is 
analyzed. In an MCG study on VT patients this 
variability is higher compared to non VT patients 
but as in TDA, QRS duration seems to be the most 
suitable parameter for risk analysis [13]. 

2.2 Intra QRS-fragmentation 

Fragmented depolarization can be detected not only 
at the terminal part of the QRS complex but 
throughout the total depolarization of the ventricle. 
Oeff et al. [14] analyzed the fragmentation of the 
QRS complex on the basis of a fragmentation index 
(FI). For this purpose the minima and maxima 
between onset and offset of the filtered QRS 
complex using a binomial low pass filter were 
determined. FI was calculated as the product of the 
number of extrema and the sum of the amplitude 
differences, increasing in VT patients. Dichotomized 
at 36, the authors found a sensitivity of 80 %, a 
specificity of 93 %, a positive predictive value of 
66 % and a negative predictive value of 96 % when 
comparing 26 healthy subjects with 32 post-MI 
patients without VT and 10 VT patients. 

2.3 Magnetic field mapping of the QRS 
integrals 

Endt et al. analyzed the spatial distribution of intra- 
QRS-fragmentation and its performance in a total of 
31 post MI patients, 20 of whom had documented 
VT using a minigrid of nine locations within the 
standard Saarinen grid. After averaging and filtering 
the signal of each channel with a digital binomial 
high-pass filter (fc = 37Hz) and a low-pass filter (fc 


= 90 Hz), the degree of curvature changes of the 
QRS-morphology was analyzed and quantified on 
the basis of two parameters: M, describing the 
number of extrema and S, describing the number of 
extrema and their heights. In general the patient 
groups could be separated using these parameters 
but there were marked differences between the 
single locations. Best values in terms of sensitivity 
and specificity were obtained for the positions above 
the middle and lower part of the sternum [15]. 

2.4 Relative smoothness score 

The relative smoothness score is used to detect 
abnormal repolarization at the beginning of the ST 
segment. The similarity of consecutive maps is 
analyzed 25 ms apart, starting at the end of QRS and 
the correlation coefficients for each subsequent pair 
of field distributions is calculated. The relative 
smoothness score represents the sum of these three 
coefficients. The lower the score the higher the 
repolarization abnormality. Until now, this approach 
was only used in small patient groups with different 
underlying cardiac diseases such as long QT 
syndrome and cardiomyopathy [16]. Data in CAD 
patients are not available, although the method 
seems promising for MI patients with risk for VT. 

2.5 QT dispersion 

QT dispersion (QTd), measured as the difference 
between the maximum and minimum QT duration in 
standard 12-lead ECG, has been proposed as a 
marker for the inhomogeneity of ventricular 
recovery times [17]. However, by taking only two 
extreme values into account, QTd is very sensitive to 
outlier error. MCG allows the simultaneous 
registration of cardiac activity at multiple sites in a 
plane above the thorax and enables the analysis of 
spatial aspects of cardiac repolarization. Fooking 
only at QTd as a temporal approach, Yukinaka could 
show that the QTd values measured in MCG were 
higher compared to those measured in 12 lead ECG 
[18]. 

We have studied temporal and spatial aspects of 
QTd in a study including 20 healthy subjects as well 
as 42 patients after MI, 11 of whom had documented 
VT (see table 2). MCG was performed using a 37 
channel biomagnetometer (Siemens Krenikon, 
Erlangen, Germany). Data were recorded for a 
minimum of sixty seconds with a bandpass of 1-200 
Hz at a sampling rate of 1 kHz. The single beats in 
the raw data of each subject were averaged and Q 
onset and T offset were determined in each channel 
to the nearest millisecond. For each subject, the 
temporal dispersion of every interval was calculated 
as the difference bet ween its minimal and maximal 



Table 2: Patient data. N = healthy subjects; EF = 
ejection fraction, see text for further abbreviation. 



N 

MI 

VT 


n = 20 

n = 31 

n = 11 

age (years) 

54 ±7 

58 ± 10 

62 ±7 

male 

6 

24 

9 

EF (%) 


61 ± 11 

52 ±7 

1/2/3 vessel disease 


17/5/4 

3/6/1 

Anterior/Inferior MI 


8/23 

9/2 


value. The spatial dispersion was examined visually 
by plotting the deviations of each channel from the 
minimal value with respect to its position (Figure 1). 



Figure 1: Three examples of the spatial distribution 
of QTd showing higher QTd values as well as a 
rougher pattern in the MI patient and moreso in the 
VT patient (see text for abbreviations). 

Healthy subjects demonstrated a very similar spatial 
pattern of QT values characterized by a line of low 
values from the upper right to the lower left. In 
contrast, the pattern of the patients showed marked 
differences in terms of local and global variability of 
QT duration. Beside higher values for QTd, local 
irregularity increased and the global pattern changed 
dramatically from that described in healthy subjects. 
For quantification we used two Smoothness Indexes, 
SI reflecting local inhomogeneity and Sin reflecting 
the overall pattern of repolarization times (described 
in detail in [19]). 

On the basis of ECG, QTd as a temporal approach 
was not able to discriminate between patient groups. 
In MCG, the patient groups could be separated but 
with a wide overlap of values (Figure 2). 

Taking spatial aspects into account, the visual 
impression could be confirmed on the basis of SI 
and Sin, demonstrating statistically significant 
higher values for VT patients. Best discrimination 
between patient groups was given by Sin as a 
marker for global distribution of repolarization. 
Using a cut off point of 5, we found a sensitivity of 
100 % and a specificity of 78 %. 


Figure 2: Values for QTd, SI and Sin in all groups 
(see text for further explanation and abbreviations) 


QTd SI Sin 



3 Discussion 

MCG offers promising approaches in the identifica¬ 
tion of post MI patients at risk for malignant ar¬ 
rhythmia. Methods such as late ventricular magnetic 
activity with both time domain and spectrotemporal 
analysis, relative smoothness index, fragmentation 
of the total QRS complex and QT dispersion follow 
the idea that the risk for arrhythmia is best reflected 
in the occurrence of inhomogeneity of ventricular 
de- and repolarization. Beside the completely non- 
invasive character of the method without any skin 
contact, the most obvious advantage of MCG com¬ 
pared to 12 lead ECG is certainly the higher number 
of registration sites. This allows the analysis of most 
of the parameters mentioned above not only under 
temporal but also spatial aspects. This increase in 
information may be of great advantage in discrimi¬ 
nating between patient groups. 

Nevertheless, a number of limitations exist. Apart 
from the small study sizes, there is no agreement 
concerning the registration technique, in part due to 
the variety of systems. On the other hand, data 
analysis must be standardized in order to allow 
comparisons and the establishment of common 
criteria such as cut-off points. This is essential in 
order to assess the true value of MCG in risk 
analysis in routine clinical practice. 
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